FOKKER-PLANCK EQUATION IN BOUNDED DOMAIN 



by Laurent CHUPIN 



Abstract. We study the existence and the uniqueness of a solution ip to the hnear Fokker-Planck 
equation — div((/3F) = / in a bounded domain of R** when F is a "confinement" vector field 

acting for instance like the inverse of the distance to the boundary. An illustration of the obtained 
results is given within the framework of fluid mechanics and polymer flows. 

Equation de Fokker-Planck dans un domaine borne 

Resume. On ctudie I'existence et I'unicite de solution ip a I'equation de Fokker-Planck lineaire 
—Aip+ div(i^ F) = / sur un domaine borne de R'' lorsque F est un champ de vecteurs "confinant" 
par exemple comme I'inverse de la distance au bord. Une illustration des resultats obtenus est 
donnee dans le cadre de la mecanique des fluides et des ecoulements de polymeres. 



1. Introduction 

In this paper we are interested in the so caUed Fokker-Planck equation 
(1.1) - A95 + div((y5 F) = /. 

• In the simplest case (that is F = 0) this equation is known as the Laplace equation (when / = 0) 
or as the Poisson equation (when / ^ 0). The solutions of these equations are important in many 
fields of science, notably the fields of electromagnetism, astronomy and fluid dynamics, because they 
describe the behavior of electric, gravitational and fluid potentials. 

• More generally, the main reason of the physical interest of equation (|l.ip comes from the fact that 
it can be put in conservative form div(J) = / with J = —Vip + ipF. Thus it can be connected 
to a generalization to the Pick's law J = — V(/? connecting diffusion flux J and concentration ip in 
inhomogeneous environments, see [9| I22j. 

• In the dynamical systems framework (see for instance |26| ) the non-stationary Fokker-Planck 
equation dt^p = eA(p — div((/9F) is usually introduced. In this case, the function (p represents the 
smooth probability density of a population driven by F and subject to e-small diffusion in the 
following sense. The term 1^9 F is a vector field representing the population ip moving with the flow 
of F, and so the divergence of this vector field represents a thinning out of the population due to F, 
which therefore contributes negatively to the local growth rate of the population, dfip. This explains 
the drive term. Meanwhile the term eAip represents e-small diffusion, and contributes positively to 
the growth rate. The study which is presented here concerns in particular the existence and the 
uniqueness of a steady-state solution. We note that the theory is closely related to applications, 
because the steady-state ip is an e-smoothing of the measure on the attractors of the flow of F 
(see |26j ) and therefore in numerical and physical experiments ip can be used to model the data 
with £-error. 

• According to the contexts, the vector field F can take various forms. In particular it may occur 
that physically realistic assumptions do not make it possible to conclude only with the already 
known results. We will give such a caricatural example in the last part. 
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Besides the problems of the existence and the uniqueness, the question which interests us is to know 
which boundary conditions are needed to ensure the existence and the uniqueness of a solution of 
equation (jl.ip in the bounded case. 

We will see that this depends on F. When F is regular enough, i.e. does not diverge too quickly 
at the boundary, data on (p at the boundary of the domain enable to ensure the uniqueness of the 
solution. We remind of this result at the beginning (in particular because the proof resembles ours) . 
We also say that when the domain does not have boundary, for instance if we are interested in the 
space R'* or on a compact variety without boundary, uniqueness is ensured by imposing the average 
of (f. 

We prove that, in the bounded case, when F is not so regular, the "good" condition to ensure 
uniqueness is still to impose the average of (p, and that in that case, the unique solution vanishes on 
the boundary. 



Except for the case where F = 0, a particularly simple case corresponds to F = WV {V is assumed to 
be regular and differentiable) and / = 0. In this case (p = exp(F) is a solution of (jl.ip in the bounded 
case as well as in the compact case or in the unbounded case. In the same way, we can easily prove 
that the case F = WV + G admits the solution ip — exp(T^) if and only if div(G) + G • WV — 0. Up 
to a renormalization, the average of such a solution will be equal to 1, as soon as exp{V) G L^{fl). 
The average value is consequently an essential ingredient to have uniqueness of the solution of 
equation (jl.ip and we could be interested in the following problem 



In the compact case without boundary E.G. Zeeman [26| proves the existence and uniqueness of a 
solution p for an arbitrary smooth vector field F (and without term source / = 0) on a compact 
manifold by the Pcrron-Frobcnius method: 

Theorem 1.1. — Let fl be a compact manifold without boundary. 
Jf F e C°°(il) then there exists a unique nan negative sohition of (|1.2p . 

Without proof in the non-compact case (as writes E.G. Zeeman on p. 152, the extension of such 
results to non-compact case is an open question), E.G. Zeeman gives some example with 17 = 
K"^. These examples show that in the unbounded case uniqueness follows from some "boundary" 
conditions on F, which are given by the behavior of F outside large sphere. Many other works 
concern these equations of the Fokkcr-Planck type in M'*. Most of these works describe specific 
assumptions for the potential F at infinity. Let us quote as an example the beautiful recent series 
of works by Herau, Nier and Helffer [13l [14] about the linear kinetic Fokker-Planck equation. See 
also the article [20] by Noarov in wich the author gives some smallness conditions in some norm and 
rapid decay at infinity for F to ensure the existence of a solution other than identical zero (in the 
case / = 0). 

In this article, we are interested in a possible generalization in bounded domains. Usually, such a 
problem is coupled with boundary conditions (Dirichlet, Neumann or mixted boundary conditions). 
For instance, the natural weak formulation of the problem 



1.1. Some known results on equation (|l.ip 



(1.2) 




Ap+ div(</jF) = / 
with p ^ 



in 17, 
on 917 
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is written 

{Find LP e Hl{VL) such that for aU V' e Hl{n) 
Jn Jn 

where (/, i/') corresponds to the duahty product i/g (il)). In order that all terms in (|1.3p 

be defined, the minimal hypotheses on data are: / G H^^{Vl) and, thanks to the classical Sobolev 
injections, F e L'^* (fi) where d^, — d for d ^ 3 and e ]2, +oo[ for d — 2. Within this framework, 
we have (see [7j): 

Theorem 1.2. — Let ^l be a bounded domain of M'', d > 2. 
If f G H~^{^) and F G L'^*{^1) then there exists a unique sohition of (II. 3p . 

The proof of the generalization which we propose is primarily based on the proof of this theorem. 
The main difficulty which appears for the study of problem (|1.3p is the following: although the 
operator A is coerciv, the operator — A + div(- F) is generally not coerciv. The reason for which the 
result is still valid lies in the (conservative) form of the term div((^ F). Let us note that an equivalent 
theorem can be proved (see [7] ) for equations of kind —Aip + G ■ V(p = f but that it is not possible 
to obtain a similar general result for equations of the type —Aip + div((^F) + G • W(p = /. In fact, 
the sum div((/7 F) + G • \/ip makes appear zeroth order terms and it is well known that the solutions 
of — Aip + A = 0, A G R, with homogeneous Dirichlet boundary condition are not unique. 
Moreover, J. Droniou proved that the same result is valid for other boundary conditions as non- 
homogeneous Dirichlet, Fourier or mixed boundary conditions (using more regularity for the domain 
n, say with Lipschitz continuous boundary). Concerning the Neumann boundary conditions, J. Dro- 
niou and J.-L. Vazquez recently showed that the same problem admits, for each fixed mean value, 
a unique solution with the said mean value (see 0). 

An other question is debated in the present paper: Which necessary and sufficient conditions must 
be placed on /, and what are the degrees of freedom on the solutions if F is not regular enough ? 



1.2. An (partial) answer 



We show in this article that if the normal component of the vector F(x) behaves |^ like 
in a neighborhood of the boundary dQ with a > 1 then there exists a unique solution ip to the 
Fokker-Planck equation p.ip as soon as the average of ip is given. Moreover, we can show that this 
unique solution automatically vanishes at the boundary 951. More precisely we prove 

Theorem 1.3. — Let ft be a bounded domain of R'^, d ^ 2. 
Let f g H^j^ (this space will be precised later) and F = k + W where k e L^{Q.) and V G C°°{Vt) 
satisfies V — ~oo on the boundary of V,. Under assumptions ( |HiD , ( IH2D and (jHaJ (see more details 
page fTS)) then there exists a unique solution of the Fokker-Planck equation (II. ip suci that Jq'P = 1. 



Obviously, the additive assumptions dTiaP enable to take into account the examples 

where the normal part of F(x) behaves like jist(x,9n) ^-"^ ^ neighborhood of the boundary dfl with 
a > 1. 

They are not satisfied when a ^ 1. Thus, there remain many cases without answers: for instance, 
when the normal part of F behaves like jist^x.an) '^i^h a ^ 1, and when the normal part of F 
behaves like dist(x^oo)'^ ^^^^ a > and /3 > 1. 

The assumptions jHil , and jHal on the potential V are rather difficult to apprehend. The 



^^"^ We can verify that we have F ^ L''*(f2), and that consequently the announced result is a generalization of the 
Theorem [L2l 
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reason for which we have these assumptions is the following: they are used in this form in each steps 
of the proof (primarily in the various lemmas) . Thus if one of the steps of the proof can be shown in 
another way that presented here, we can hope to be freed from certain corresponding assumptions. 



1.3. Outline of the paper 

The paper is organized as follows: 

• In Section O we give mains tools adapted to the studied problem. First of all, some tools about 
differential geometry to understand "explosive" boundary conditions. Next we give all the lemmas 
which are used in the main proof. 

• In Section [3] the precise statement of the main result is enunciated, see Theorem 13.31 page [TBI 

• The Section [5] is devoted to the proof of Theorem 13.31 It is composed of two parts: the existence 
proof and the uniqueness proof. 

• The last section (Section [5]) gives an application to fluid mechanics and some numerical results. 

2. Main implements 

In [5], the author gives an existence and a uniqueness result for a Fokker-Planck equation for a 
particular vector field F and in a particular domain Q which is a ball. For more complex domains, we 
must understand the effect of the geometry in the proof. We will present in this part some elements 
of differential geometry adapted for our calculus. Next, we will precise the functional framework 
adapted to the Fokker-Planck equation of this paper. Finally, we will give multiple fundamental 
lemmas which are use for the proof of the Theorem 13.31 



2.1. Elementary differential geometry 

The results of this part are largely inspired on Subsection 2.1 of the paper [2] and on the annexe C 
of the book [3] . 

Let n he a smooth (say C^) bounded domain in W', d ^ 2. We denote by F its boundary and by v 
the outward unitary normal to F. The distance between any x G M.^ and the boundary F is denoted 
by (5r(x). 

For any e ^ 0, we introduce the open subset of fl: 

Oe = {x e f2 ; (5r(x) < e} and fi^ = {x e f2 ; ^r(x) > e}. 

It is classical that, for e small enough, the two maps Sr (called distance to F) and Pp (called 
projection on F) exist and are regular on O^- This allow to use tangential and normal variable 
near F, defining for any function / : M the corresponding function / : [0, e] x F ^ K bj{^ 

y{r,9) e [0,e] X F, /(r, 0) = f{9-rv{e)). 

Moreover, for any / G L^{0^) we have the following change of variables formula: 

(2.1) / /(x)dx= / r f{r,9)\J{r,9)\drd9, 

Jo, JrJo 

(■^^ For non-ambiguous cases, / and / will be together denote /. 
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where dr and d6 corresponds to the Lebesgue measure on [0, e] and F respectively, and where J 
is the jacobian determinant of the previous change of variables. Introduce km — supggp supj Ki{0), 
where ki{6),..., k,(1_i{6) are the principal curvatures of T at 0, we have for ^ r ^ e and for S F 

{l-ennf'^ J{r,0) ^ 1. 
Using the change of variables (|2.ip we deduce that for / e L^{0^) we have 

(2.2) (i-sKMf-^f I f{r,e)drde^ / /(x) dx / /'/(r, 0) drd0. 

JrJo Jo, JrJQ 

Roughly speaking these relations show us that for thin tubular neighborhood of F in the ja- 
cobian J(r, 6) is equivalent to 1 for small r, uniformly with respect to G F. Consequently, the 
relations (|2.2p give an approximation of / by /p f. 

Notice too that it is possible to define 5t as a regular function on ^l. This extension will be also 
noted 6r. 



2.2. Functional spaces 



The Fokker-Planck equation (|l.ip make appear a potential V, assume to be smooth on il, via the 
force F = K + WV. This potential is supposed to be confinant so that we assume that G L^{fl). 
We can always define a maxwellian function M by 

Notice that M g C°°{^1 : M*^_) with J^.^ M = 1. Notice too that the interesting cases for the present 
study correspond to a maxwelian M vanishing on the boundary F = dfl (that is F = — oo on F). All 
the maxwellians considered in this paper will satisfy M ~ on F. The Fokker-Planck equation p.ip 
is written 

(2.3) div(^/^-Mv(^)) = /. 

From the peculiar form of this Fokker-Planck equation (j2.3p . the adapted functional spaces use 
Sobolev weight spaces on the domain C M''. More precisely, we introduce 

LIj ML^{Mdx) and H\[ := MH^{Mdx), 

endowed with their natural norms respectively given by 



/ M 




2 


In 


M 




1 M 




2 

+ M 


Jn 


M 





Note that a large literature on weighted Sobolev spaces exists, see for instance [231 US], the references 
given in the notes at the end of Chapter 1 of [12] or the classical book [51] . Among all the weights 
which are generally considered in the literature, only some are vanished on as it is the case 
of the weight M . We refer to [21] for references of spaces with weights which equal to a power of 
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the distance to the boundary. Nevertheless, the spaces L'^{Mdx) and H^{Mdx) being spaces with 
"traditional" weights, the spaces and H^^ are not it. For this reasons, the results shown in 
the following section are not always "traditional" corollaries (except some whose proof is a direct 
consequence of the results in L'^{Mdx) and H^{Mdx), see for example the proof of lemma [2T4|) . 
Note first that the two spaces and H^j are Hilbert spaces (see [Ml Th. 3.2.2a]). We introduce 
the following normalized subspace 

Hli,o ■■= &Hlr, / ^ = 0}. 

Jn 

In the sequel, the space ff^ will by equipped too with the semi-norm defined by 




and we will see (Lemma 12.91 on page [T3l) that || • Hhj^^^ is a norm on the space H\j q. Finally, we 
denote H^j the topological dual of H\j g, that is the set of continuous linear forms on H\j q. Each 
application x G ^a/ ^'^^"^ defined by x : </? G H\[ q ^ (x, ^) G M. By its continuity, for each 
X G H'^j there exists C G M such that 

As it is usual, the smallest of these constants C is denoted ||xll^f-i- it is the norm of x on H~^^ . 



2.3. Properties of the functional spaces 



The proof of the main theorem (Theorem 13. 3p follows the ideas of J. Droniou Nevertheless, the 
proof of J. Droniou, given in the case where F G i'^*(r2), does not use these "degenerated" spaces 
and use traditional results concerning usual Sobolev spaces. The essential contributions which are 
presented ties in the fact that these "classical" lemmas in the case where M does not vanished on Q. 
are still true (sometimes in a weaker form) when M is a maxwellian as previously introduced, and 
in particular when M = on F. So, the goal of this subsection is to give some essential properties of 
these functional spaces (Poincare-type inequality, Sobolev injection, compacity result. Hardy-type 
inequality...). 

Notice that in the estimates, the symbol < means "up to a harmless multiplicative constant" , allowed 
to depend on the domain only. 

The first result that we present is a result allowing to controlled in {Q) as soon as ip £ H\j . 

This result can be seen as an inequality of Hardy- type We will prove this first lemma under the 
following assumption!^ 

VflM(O) = 0, 



(Hi) 



Op 



(•^^ Hardy inequality indicates that for a function / defined on a bounded domain Q of and vanishing on the 
boundary F of £7 we get \ f I &v\h'^{Q:) $ I^/Il^jq) where (5p corresponds to the distance to the boundary F. 

Recall that, as it is specified just before, in this article the function M is a Maxwellian function which vanished 
on the boundary of the domain. The assumptions introduced here are consequently additional assumptions. 
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These hypotheses will be only used in the neighbourhood of the boundary F of fJ. In such neigh- 
bourhood, the notation V/j corresponds to the radial derivative (that is to say the derivative in the 
direction of the normal vector to the boundary) . 

The various assumptions introduced in this part will be discussed in Section [3l Nevertheless, it is 
important to notice that the three assumptions formulated in jTlil ^-re independent. For example, 
in the radial case the function M defined by M{r) = y^r satisfy the last point of Pi.il but does not 
satisfy the second point. Reciprocally, M(r) = r/ln{r) satisfy the second point but does not satisfy 
the last point. 

Lemma 2.1 (Hardy- type inequality). — If M satisfies \Tii\ then for any Lp e iJ^^ we Jiave 







In "r 


M 



< 



In fact, as can be seen from the proof, we shall not need entire norm of ip in H^j since we will only 
use the radial part of the gradient. 

Proof — For any e small enough (more precisly, e ^ en, see Part 12. 1|) we have ft — H.^ U Oe and 



M 


V 


2 r M 




2 f M 




2 


51 


M 




M 


+ / 6^ 


M 





Since l/dr is bounded in we easily deduce that 



r M ^2 ^ r 



M 



M 



and the main difficulty to prove the Lemma 12.11 is concentrated in the control of the integral 

2 



r M 






M 



For similar reasons, we can suppose that 
(2.4) 



0. 



In fact, introduce a regular function 7^ (only depending on 17 and e) wich vanished outside the 
e-neighborhood of F and is equal to 1 in the |-neighborhood. The estimate in Lemma [??T] clearly 
holds for (1 — 7e)iy9 and the proof therefore has to by conducted only for 'ysf-, which for sake of 
simplicity will be denote p in the following. 




The value of e being able to be chosen as small as desired, we shall use copiously the change of 
variables (|2.ip as well as approximation (|2.2p . 
For each e F let introduce the function 



Note that the relation (|2.4p implies that /ie(e) — 0. 

To simplify the notations, when there are no ambiguities, we do not note the dependence with 
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respect to the variable in he and in Ad . The integral /q writes using the approximation 



where s = 



'r Jo ' 

The goal is to control Iq with ||<y3|j/fi^. We will proceed in two steps: 

(1) We prove that h{r)M' {r)/M{r) = on r = 0; 

(2) We control Iq with ||iy9||^i^. 

• Step (1): We use the following change of variable adapted to the maxwellian M: 

$ : (r, e) e ]0, e] X r I — > {s, 6) e [0, +oo[xr, 

dt 

M{t, e) ' 

Notice that the jacobian determinant of $(r, 0) egals — 1/M(r, 0) and therefore, M being positive 
on n, it does not cancel on ]0, e] x F. Moreover, for all G F we have limj.^o ^{f, &) = (+oo, 6) since 
using the assumption on M we get dr/M{r, 6) = +oo. Consequently $ is a local diffeomorphism 

from ]0,e] x F to $(]0,e] x F) = [0,+oo[xr. 

For any function / : ]0, £[xF i— > M we will define the function, noted / : ]0, +(X)[xF, by / = / o $. 
Using the change of variables introduce in Subsection 12. li for any function / : i— > R we have 
define the functions / :]0, £[xF i-^ R and / :]0, +oo[xF such that, with the previous notations for 
the name of variables: 

/(x)=/(r,0) = /(s,0). 
In the nonambiguous cases, we will note /, all the functions /, / or /. 

As announced before the proof, we only need the radial part of the gradient in the desired estimate. 
In term of new coordinates, the radial part of the gradient of a function / defined on fi^ corresponds 
to the derivative with respect to the variable r in the new coordinates: Vij/(x) = drf{r, 0). 

For any g F let ge be the function defined on ]0, +oo[ by 

Derivating with respect to the variable s, we obtain (as previously the variable 9 will be understood 
as a parameter and we do not note its dependence): 

(^(x) ^ ^ 



M(x) 



dx^-\g'{s)\^\J{r, 6>)|d6» ds. 



We deduce that (using the approximation J ^ 1 vahd for e small enough) 



M(x) 



l[l de. 



VM(x), 

Since ip G H\,j the integral Ii is finite. Consequently, for almost every G F we have 

p+co 

I \g'{s)\'^ ds < +00. 
Jo 

We deduce that 

Va>0 3sa>0 / |5'(s)p(is ^ 
Next, using the Cauchy-Schwarz inequality we obtain, for any s G ]0, H-cxd[, 



7(s)=g(s„)+ / g'it)dt^ 



g{sa) + ay/s — Sa if s > Sq, 

sup g if < s ^ Sa- 
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We deduce that for almost every e F, for all a > there exists a constant Ce.Q > such that, 
for any s G ]0, +oo[ we have g{s) ^ a^/s + Ce.a- Using the r variable, this result is written: for any 
r G]0,e] 



Arjr) 
M{r) 



h{rY a^M'{r) 



dt 



that enables, see assumption to obtain for almost every e F the relation 



(2.5) 



M'{r) 
M{r) 



h{rf 



0. 



r=0 



Step (2): Now, we prove the lemma. Since Lp G i?^/, we know that 



h 



Af(x) 



( y'(x) 



dx < ||</5||^l < +00. 



^Vm(x) 

We express Ii making appear the h function and using the change of variables x i— > (r, 0) together 
the usual approximation for the jacobian determinant of this change of variable (see the Part 12.11 
and the relations (j2.2p ). We obtain 



h = 



where, for e F, the quantity li{0) is defined by 



(h\rf 

Moreover, an integration by part gives 

Jo Mir) ' ^ ' 2 

We deduce that 



M'ir) ^, 



2M{r)./M{r) 



M(r) 



h'{r)h{r) 



AM{ry 



h(r 



dr. 



M'{r) 
M{r) 

h{rf 



h{rf 



I {h\rr + X{r)^) dr 



-M'(r) 
. M{r) 

h{r? 



with 

The assumption ( |?iiD on M is written A — a > — 1. Moreover using the equation (|2.5p and the 
relation (|2.4p the braket term is vanished. We obtain 

h{rf^ 



((M'{rW 




\\ M{r) ) 





'h{r)^dr^^ 

M'(r) 
M{r) 

2 



hi0)>J^ [h'ir)'-a- 



4^2 



dr. 



Moreover, thanks to the Hardy inequality (holds since h vanishes at 0, this is a direct consequence 
of equation p.Sp ). we deduce 



(2.6) 



Ii{e) ^ (1 - a) / h'{rfdr. 
Jo 



Since a < 1, this control allows us to estimate h'{r)^dr. From the Hardy inequality again, we 
obtain the following estimate 

^2 



r ^^dr 4 h'irfdr ^ (0) 
Jo Jo 1-a 



Integrate with respect to the variable 0, we obtain 

h{rf 



lo 



dr de sC 



r Jo 



1 -a 
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r M 


f 




M 



The function Lp being in H\j^ right hand side termes are bounded by ||t/^||?Ti • Up to the change of 

M 

variables we deduce 

which imphes the announced result. □ 
The additional hypothesis 

(H2) 3c> |VM| 

c Or 

implies from Lemma [^?T] that ii (p E then '^^^= e It is in this form that the Lemma [^TT] 

will be generally used in this article. 



uiuj l|^.U|l aiiuw UD Liiau vj^y-^^ uciuiigo uu j_ ^ 

propertie is completed by the next lemma: 



Note that in term of x variable, the inequality (|2.6p show us that V belongs to This 



Lemma 2.2 (Inclusion). — If M satisGes IjfHil and ( [7^2 D t-hen we have the foUowing inclusions 

LliCL^iQ) and C H^in). 

More precisely, if Lp £ L\j then pj^fu G and if p e H\i then pj^fu e H^{n). 

Proof — Inclusion C L'^( fl) is obvious since on the one hand, by definition of L\,j, we have 
(/? e if and only if ^p/VM G and on the other hand \/M G L°°{n). 

To prove the inclusion Hlj C 11^(^1) we use the Lemma [2.11 with additional assumption jH-zl '- if 
(fi G then we have 

1 VM ^ ^ 



Va/M/ VM / VAf/ 2 M ./m 




Consequently p/y/M G H^{fl). Hence this function ip/^/M has a trace on the boundary F. Since 
M is a regular function on fl vanishing on F, we deduce that p = -^=\[M G □ 

This next lemma is interesting in themselves for understanding the space better. Moreover, it 
will be used in the proof of the Lemma [^771 



Lemma 2.3 (Density). — If M satisfies ( |?iiD and then we have the following equality 

ttI 

Proof. — Let p) G H\,j and define, for n G N*, the function pn by 

{1 ifOs^<<l, 
2 - t if 1 sC < < 2, 
if i > 2. 

We successively prove that 

(1) the functions Pn are in H]^j, 

(2) we can approach these functions pn with functions, 

(3) the sequence {</3„}„gN. converges to p) in H\j sense. 

These three points clearly implicate the lemma. 

• (1) By definition of x, we have for all n G N* the relation \pn\ ^ \p\ hence ^f" |^ ^ /f2 -'^'^ I "H" I 
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which is bounded since Lp E H\j C L\j- 
we write 



To control the gradient part of the iJ^^-norm of iy9„ , n G N 



i 



Mn\nM) ' M 





+ 







I 1 y VA/ 
I n M M2 



M \^\nM 
term is 

and can be controUed by 2 



By definition of the trucature function x, the last term is not egal to if and only if 5^ ^ < 
In this case it write 



M M I 





2 




2 


f 4 


V3 VM 


VM/ 




\m) 






M M 



We have 
2 



Since G i? and using the Lemma 12.11 together with the assumption (jHg]) , we deduce that 
/j^M|V(^)| is bounded. Consequently, (p„ G H\j. 

• (2) For each rt e N* the function (y3„ e C iJ^ is egal to in a neighborhood of the boundary 
r. Approaching (/3„ by a sequence m}mGN (which is egal to on a neighborhood of V) in 
allows to approach ipn by the same sequence {fn,m\mefi in H\j. 

• (3) Note that (y5(x) — (p„(x) egals to if A/(x) > i. In the other case, that is for all x e il such 
that M(x) ^ i, we get |v5n(x) — V3(x)| < |</5(x)| since < x ^ 1- Hence 

2 



/ M 




\ f M 


V 


In 


M 


Jxen;M(x)^i 


M 



Since (/? € C i^/ and since the measure of the set {x G ; M (x) ^ i} tends to when n tends 
to +00, we have 



(2.7) 



M 



M 



n — *+oo 



0. 



Concerning the gradient part, we write 



v(x 



(tim) 



The first term is non zero if Af ^ i and is bounded by | V(-p^) | in this case. The last term is non 
zero if ^ ^ M < i and is bounded by 2 1 I . Hence we have 



M 



M 



-0 


2 




VM 










Af 


/M 



/xGO; Af(x)^-l- 

Since ip G H\j, the first term is bounded, and using the result of the Lemma 12.11 again, we know 
that the last term is bounded too. More precisely, as previously, using the fact that the measure of 
the set {x G ; M(x) ^ tends to when n tends to +00, we have 

(2.8) 



M 



M 



■>+oo 



0. 



The relations (|2.7[) and (|2.8[) ensure that the sequence {fn}n£n' converges to ip in H\j. □ 
Now we introduce a compacity result for the spaces and H\^j which is comparable to the classical 
compact injection Hl{Vl) ^ L^{il). 



Lemma 2.4 (Compacity). — The injection Hi 



M 



L\j is compact. 



Proof. — To prove this lemma, we use the following result due to G. Metivier [HI Proposition 3.1 
p. 221] affirming that the weight Sobolev space injection H^{Mdx) C L'^{Mdx) as soon as we 
have M>0onri,M = 0onr and VM ^ out of F. Notice that this last point is a consequence 
of the two first points M > on and M = on F, at least in a neighborooh of the boundary F, 
what is sufficient for our results. For sake of simplicity, we will nevertheless work with fl. 
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Consider a sequence {<Pn}neN bounded in H^j and show that a convergent sub-sequence can be 
extracted. By definition of H^j, for aU n e N, there exists gn € H^{Mdx) such that (pn = Mgn- The 
sequence {iPn}n£N being bounded in iJ^^, the sequence {gn}n£¥i is bounded in L^. We use here 
the result of G. Metivier. We can extract from the sequence {gnjnGN a sub-sequence, still noted 
{5n}nGN ^nd such that 

gn ^ g in H^{Mdx) and gn ^ g in L^{Mdx). 

By definition of the spaces H^j and L'j^,J we conclude that 

ipn Mg in and <^„ ^ in , 

which proves that the injection H\j > L^j is compact. □ 

In the same way, we present the next lemma which proves that functions in H\j are in certain L\j, 
p > 2, where the weighted-space is defined by 

p\i/p ) 



Ll, = {^eLUn); 



M 

n 



M 



< 



and endowed with its usual norm. 

This kind of result is essential for the proof of the main theorem (Theorem 13. 3p ; it is proved under 
the assumptions and (jHg])- We will note that assumption is used in the following weak 
formulation (obtained by integration): 

m^i 3c' > M ^ c" sy^ 



More exactly, we have 



Lemma 2.5 (Sobolev-type injection). — If M satisfies ( |HiD and ( [7^2 1 ) t-hen there exists p > 2 
such that the injection H\j ^ is continuous. 

Proof. — First, let us note that ip e L^^ if and only if ^^^^^^ G L^(17) where the spaces LP{il) 
are the classical Sobolev spaces on the set fl. Let ip E Hm- In the next three steps we will prove 
that there exists p > 2 such that — -. — rrr- G L^ifl). 

• Since p G we clearly have tp £ L\j and so pi/^/M e L^. Moreover, using assumptions 
Lemma |2. II with assumption we obtain 



ip \ _ V(p VM (p 2 




e L'{n). 



Consequently, ipj^fM E iJ^(ri) and using the classical Sobolev injections, we deduce that ip/VM E 

iy'(r2) for all q ^ 2d/ {d — 2) (and for q < +00 in the 2-dimensional case). 

• Using the assumption (jWg]) and the Hardy inequality (Lemma 12. ip we obtain 



VMM'' Sr Vm M^ 
• From the two previous steps, we can write that for any /3 G R such that ^ /? ^ c we obtain 
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with r = q^^2(c-;3) • note p the real such that 1 ^ ^ = /i^ + ^- The previous result is written: 

for any p G M such that 2 ^ p ^ T^2c obtain 

^ eL^ with ^^^'^ 



iV/i-i/p qp - 2q + Acp - 2p + i' 

In particular, we have Zi/p ^ 

as soon as r ^ p. The inequality r ^ p holds if and only if 
p ^ 2 + 1^^?^. It is thus possible to find p>2 such that — ^^-^ e i^'. □ 

Remark 2.6. — According to the previous proof, we obtain the inclusion HIj C for all 
p ^ 2 + ^^pzfjqij- For instance, in the two dimensional case {d = 2) the inclusion C holds 
for all p ^ 2 + 4c. 

To build functions in H^j g we use the following lemma which will be important to obtain a lot of 
test functions in the weak formulation later. The hypothesis jTii^j allows to use the Lemma [2731 



Lemma 2.7. — Assume p-C2l hold. Let ip e HIj and ^ : M ^ M be a continuous application, 

piecewise-C^ such that ^' is bounded on R. Then we have 

„-,iv(|^)=f(|-)v(^)a„d|b||„.,„S||f|U-,.,||*||„;„, 

Proof. — The proof of this lemma uses the Stampacchia lemma which affirms that if g G H^{uj), 
io being an open subset of , and ^ : R ^ R is continuous, piecewise-C^, such that ^' is bounded 
on R then we have i{g) S iJ^(w) and V^(5) = C{g)Vg. 

The Stampacchia lemma is a local result, hence applied to g = ijj/M which is in Hl^^{Vl) it shows 
that the formula V^(-^) = '?'(m)^(m) t'^^l^s in From this formula, it is obvious that 

If G H]^. The fact that (p is null average is then immediate since j^M — 1. □ 
One more important ingredient in our study is the following linear operator 



/:^ = -div(Afv(^) 



on the space L\j and with domain, see [181 Remark 3.8, p. 9] given by 

D{C) = {,peHl, ■ ^-l|div(Mv(^))|'<+oo}. 

We also find in 18, Proposition 3.6, p. 8] the following result and its proof which will be used to 
introduce the Galerkin approximation method later. 

Lemma 2.8. — TJie operator C is self-adjoint and positive. Moreover, it has a discrete spectrum 
formed by a sequence {(n)neN such that in tends to +oo when n tends to +oo. 

Concerning the uniqueness results for a linear operator, it is known that the eigenvalue 0, that is 
the kernel of the operator £, is particularly important. 

Lemma 2.9. — The kernel of the operator C is the set {XM, A £ M}. 

Proof. — This lemma is an immediate consequence of the following formulation of the operator C: 
(2.9) 

where (•, •)^2 corresponds to the scalar product subordinated to the norm || • ||^2 on L'^.j. In fact, 
let V' be a function such that £-0 = 0. We obtain (£-0,-0)^2^ — and the formulation (|2.9p yields 
V(V'/M) = 0. Thus, thanks to the connexity of fi, we deduce that ip — XM with A G R. □ 
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The last lemma is a generalized Poincare inequality adapted to the weighted spaces introduced 
before. To obtain such a lemma, we use the fact that the potential V = ln(M) is concave. More 
precisely we will suppose that 

(H3) 37 >0 v(^) s=:-7ld. 

Lemma 2.10 (Poincare- type inequality). — If M satisGes (jHs) then for all ip G H^j we get the 
following Poincare-type inequality 



1/m 

7 Jn 



For the free- average functions (that is for ijj E Hm q) this Lemma [2.101 show that the two norms 
II • 11^1^ and II • ||_ffJj-Q this space are equivalents. This equivalence will be usually useful in the 
remainder of the paper. 

Proof. — Let ip G H^j and introduce the non-stationary problem 

( utit,x) + Cu{t,:x.) =0 for (i, x) e R X fi, 
1 ?i,(0,x) = </'(x) for X e fi. 

The following time-dependant functions 

u 



M 



M 



satisfy 



D\t) - 2 / —-ut 



and Hi 



M 







/ M 




In 









/ M 




In 





Moreover we have 



H'it) 



f-V 



For clearify the following computations, let us introduce the duality operator £* such that 



1 

M 



div(AfVw) = -Aw 



VA/ 



We have £u = MC*{^) and 



H'{t) 



MV 



(s)-(^-(s) 



Since V{C*v) = C*{Vv) - ■ Vv we deduce that 



H'{t) 



y(™^ 

V M 



) 



\m) 



The first term of the right hand side is written —2 /q V^(-p^) |^ and is non-positive. Using the 
assumption ( [TiaD , the last term of the right hand side is controled by —27 /q -^|V(-p-) |^, that is by 
-2-fH{t). We obtain 

H'it) < -27 H{t). 

Hence Hit) ^ H{Q)e~'^'^* . Integrate in time, we obtain: 

.+00 ^ 

(2.10) DiO) - Di+00) = 2 H{t) dt -HiO). 

Jo 1 

To evaluate Di+00), we consider a stationary solution u^o- We note that due to the spectral prop- 
erties of the operator C (see Lemma l2.8p . for any initial data, u tends to a stationary solution Uoo as 
t +CX3. By definition it is in the kernel of C and following the Lcmma l^?^ there exists a constant A 
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such that — AM. But the evolution equation on u implies that the mean value u is conserved: 
""cc — Jfi ^, that allows to obtain the constant A = J^^- We deduce that 

Consequently, the inequality (|2.10|) corresponds to the following one 




which exactly is the inequality announced by the Lemma [2. 101 □ 
Notice that it is possible to obtain a proof of this Poincare-type inequality by contradiction, see 
for instance [TBJ p. 7], or peraphs using the hole-space case (for example for f2 = R'^) proved in 
H.J. Brascamp [4] (see also Proposition 2.1 in [6]). 



3. Statement of the main theorem 

3.1. Definition of weak solution 

When we consider the Fokker-Planck equation p.ip with vector field F decomposed as the sum 
F = K + WV where k G L°°{n) and G L^{fl), we can introduce the maxwellian function M 
by M — -j—^ and rewrite (jl.ip as div((pK — AfV(-^)) = /. If we look for a solution with given 

average, that is for instance a solution such that Iq^^ — 1, then we can reduce to the case where Lp 
is free-average exchanging into Lp — M and / into / — div(M n). We obtain the following problem 

div(¥.K-Mv(^)) = / inf!, 

with / tp ^ Q. 
Jn 

Using the adapted spaces introduce in the previous part, the weak formulation of this equation is 
written: find p G H^j q such that for all ip G H^j q 



where (•, •) denote the duality brackets between i?^/ and HIjq. 



3.2. Assumptions on the potential 



In this article we are interested in the case where the vector fields F quickly explodes near to the 
boundary. The fact that F is decomposed as a sum of two terms makes it possible to describe all 
the "explosive" behavior in the part WV. In addition to the fact that V equals —oo on P to ensure 
the explosion, the assumptions given on V (or on M, which is equivalent) can be checked only in 
a neigborhood of the boundary P. More precisly in order to use the lemmas proved we will use the 
following assumptions 

2 /VkA/n 



(Hi) 



M 

VflM(O) = 0, 
3& > VrM 



51 



1 

M 



< 6, 
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(H2) 3c>0 |VM|s;-— , 

c Or 

(H3) 37 >0 v(^) s=:-7ld, 

where we recall that Vij corresponds to the normal derivative and where Sr represents the distance 

to r. 

Notice that we can rewrite these assumptions in term of the potential V (wich is given with respect 
to the maxwellian M by F = In M) , see for instance Theorem 13.31 page [1^1 It is important to note 
that these assumptions are satisfied for the radial functions M (i.e. functions depending only on the 
distance to the boundary) on the following form near to the boundary 

M(r)=r" with a > 1. 

In other words, the result is shown for vector fields F whose the normal component explodes like 
with a > 1. 

Remark 3.1. — As it was announced as introduction, an interesting case corresponds to the 
following Fokker-Planck equation 

-eAtp + div(</jF) /, 

making appear a small parameter e. We can come back to the previous case using F = ^F. We note 
that if we define a Maxwellian M such that F = WM/M then the Maxwellian M adapted to F, 
i.e. such that F = '^M/M, satisfies M = C M^l"^ . The assumptions on M can thus be interpreted 
on M and we show that they are less constraining in the following sense: they are checked when the 
normal component of F behave like ^ for all a > e. 

Concerning the assumption on the "interior" part /t of F = k + VV^, that is about k E L°"{n), 
we can note that this assumption is stronger than that announced by J. Droniou in [7 . In fact, we 
will see during the proof that the regularity required on k comes from a product lemma. Roughly 
speaking, if the product of a function H^{n) by a function L^'(ri) is a function L^(ri) then the 
theorem is true as soon as k belongs to LP{il). In the classical case the usual Sobolev injections 
H^{Q) C L^'*/'^'*^^) imply that p = is sufficient. In our case the injections of "Sobolev" type (see 
the Lemma [2. 5 1) are not also "generous" and a product HIj x LP{n) will not belong to L^(ri) for 
as many values of p. We can possibly improve the result of the theorem by taking k € U'{^) with 
p ^ d+l/c. 



3.3. Main theorem 

We prove in Part |4] the following theorem. 

Theorem 3.2. — Let fl be a bounded domain of M'', d^ 2. We denote by T its boundary which 
is assumed to be of class C^. Let f e i/^/ and F = k + W where k e L°°{n) and V G C°°(fi) 
satisGes V = —00 on T. 

If the assumptions jHil , ( 17^2 1 ) and jHsl hold then the problem (|3.1|1 admits a unique solution 

LP e Hln^. 

We can deduce - see the link between a free-average solution and a solution with given average on 
Subsection 13. II - the following theorem where we recall all the assumptions 

Theorem 3.3. — Let fl be a bounded domain ofM.^,d^2. We denote by T its boundary which 
is assumed to be of class C^. Let f e H^/ and F = n + \/V where k e L°°{n) and V e C°°{n) 
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satisfies V ~ —oo on T. 

If we assume that, in a neigborliood of the boundary T, we have 



(Hi) 



3a <1 (Vi^y) + 
VjiV = on T, 



3b >0 VrVc^ 



<b. 



3c> 



ivyi 5$ 



5v' 



(Hs) 37 > VV s; -7 Id, 

where Vjj corresponds to the normal derivative and where Sr represents the distance to T, then 
there exists a unique (weak) solution of the Fokker-Planck equation 

-Atp + div{ipF) = f in fl, 

such that = 1 • 



4. Proof of the Theorem [Q] 
4.1. Existence proof in Theorem 



Principle for the existence proof of Theorem 13.21 - The maxwellian M satisfying the assump- 
tions dT^iD , <Pi.2l and ( IH3D , we use the different lemmas proved in Part[2l For instance, using the 
equivalence between the norms || • and || • W^^^ ^ on the space 11]^^^, see Lemma [2.101 the operator 
— div(MV(-p-)) is coerciv on i?^Q thus we can (see for instance the Lax-Milgram theorem) prove 
that there exists a weak solution (that is belonging to Hlf q) to equations like 

^a„(Mv(^))./ 

as soon as the source term / belongs in H^j^ . Moreover in this case we have ^ ll/ll/j-i- 

Because of the non-coercivity of the operator — div(MV( -jg)) + div(- k), we start by studying an ap- 
proach problem. For each n e N, let us consider the application r„ : r e M i-^- max(min(r, n), — n) G 
R and let us denote by Fn the following application: F„ : -0 G t-^ ip E q C where is 
the weak solution of 

(4.1) - div(Mv(|:)) =/ - div(AfT„(^ 



For -0 e LIj we have MTn{ip/M) e L|j and since t!U k € L°°{Vl) we get div(Afr„(;^)/t) G H^l . 
The function F„ is then well defined. 

Let us prove that Fn is a compact application by showing that its image Fn(L\i) is bounded 
in H\j. Consider — Fn{4') € Fn{L\j). Taking i/; as a test function in the weak formulation of the 
equation (|4.ip we obtain 



M 



a 



Here the assumption on k. is essential. Following the proof of J. Droniou 7 it is possible to improve this assumption 
using the Sobolev injection 12. 51 more finely. See discussion concerning the assumptions on re page I16l 
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In other words, by using the duahty definition and the Cauchy-Schwarz inequahty, we have 

ll^llffi „ ^ll^ll//l,„ + ll'^IU~(a)i 



M 



(I) 



M 



point, denoted by ipm in L\j. This fixed point is consequently a solution of 



Using the fact successively that for all r e M we have |T„(r)| ^ n and that Jq^M — 1 we deduce 
that 

llV^llffi^^^ < U\\Hl,„+nM\L^(n)\mHl,,- 

Consequently we have 

Thus, the image of L\j by the application Fn is contained in the ball of H\,j of radius l + n||/«||ioo(Q) 
(up to a multiplicative constant depending on 17, which appears in the symbol <). Moreover, the 
injection H\j ^ L^j is compact (see Lemma 12. 4p and the application Fn is clearly continuous. 

Applying the Schauder fixed point theorem, we conclude that the application Fn admits a fixed 

-2 

■'M- 

.M - 

for all test functions tp g q . 

The continuation of the proof consists of obtaining estimates on these functions V-'n in order to be 
able to pass to the limit when n tends to +oo. 

Estimate of M ln(l + |i/>„/Af |) in H]^ g-norm - Let ^ be the application from R to R defined 
by £,{t) — Jq (i+\s\y^ ■ This application is continuous, piecewise-C^ and with a bounded derivative. 
According to Lemma [2.71 we can choose ip = M^(V'n/Af) — M J^M£^{ipn/M) as a test function in 
formulation 

• The first of the three terms obtained is treated in the following way 

2 



(4.2) 



MV 



MTr 



(l7)''-(l^)^</-> 



(4.3) 



MV 



(til 



V 



(t) 



M ^ — 



M In (l + 




) 




M 





• For the second term we obtain 



MT„(|).^V(« 



MTn{^) V(tr) 
— r-K ■ , , 



L«=(0) 



n 1 + 

M 



1+1^ 



M , 



^ \ M \ 



V(ln 1 



M 
Ipn 



M 



Using the fact that for all r e M, we have |T'„(r)| ^ |r|, we deduce thal|^ 



(4.4) 



< 



Mini 



M 



Hi 



We also use the Cauchy-Schwarz inequahty to show that 

Mf^J I M, 



Mp 



Mp 
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• For the last term, using / e -ffj;/, we deduce 



< 



' ^ "'-'-MX) 



(4.5) 



M- 



v(tf) 



'\ + 1^1 



The three estimates (14.31) , (14.41) and (14.51) enable us to obtain for all n e N 



(4.6) 



M In 1 



M 



< 1. 



Hi 



Estimate of /.t({Q E ^2 ; |i/'n(Q)| ^ ^-^(Q)}) - In this paragraph, we control the size of 
the set where ipn has a large value, that is the set = {Q G ; |'0n(Q)| ^ fcAf(Q)} for fc e N. 
The natural measure in the present context is the measure d/i = M(Q)dQ (dQ being the classical 
Lebesgue measure on il C M''), which enables to take into account the weight of the Maxwellian M. 
Writing £fc = {Q e 17 ; (ln(l + |V'n(Q)/A^(Q)|))^ > (ln(l + k))^} we obtain 



We easily deduce the following estimate 



Af ln(l 



M 



))% / Af(ln(l + |^|))% / Af(ln(l + fc)) 



Introducing the measure ~ Af (Q)dQ this inequahty is also rewritten 



1 



M In 1 



M 



(ln(l + fc))2 

Taking into account the estimate (14.61) . the previous equation is written 

1 



(4.7) 



A*({Q e ^ ; IV'n(Q)| ^ fcA/(Q)}) 



< 



(ln(l + fc))2- 

Estimate of MSk{ipn/M) in ifj^ g-norm - Recall that for fc e N the application Tk is given 
by Tfc : r e K max(min(r, fc), —k) G M. We now define the application 5*^ such that Tk + Sk = id. 
To obtain an estimate on ipn we successively obtain an estimate on MSk{ipn/M) and then on 
MTkiipn/M) for a sufficiently large keN. 

Let A; e N. Taking ip = MSkiipn/M) - M MSk{ipn/M) as a test function test in ((42| . According 
to Lemma 12.71 this choice is possible and we obtain 

"-(l)--(*(l))-//-»(l)--(^'(l7))^<^->' 



C 



A B 

• Since Sk + = id and for all r e K we have 5^ (r) = or (r) = we deduce that the first 
term A is written 



(4.8) 



A 



M 



V 5, 



MSki — 
M 



Hi 



• Using the fact that for all r e R we have |T„(r)| ^ |r| and using the Cauchy-Schwarz inequality, 
we estimate the second term B in the following way 



\B\ < \\k\ 



M 



M 
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However |7/.„/M| = |Tfc(V'„/M)+5fe(^„/Af)| s: fc+|5fc(V„/Af)| thus |V„/VM| < k^+^I\Sk{il^n/M)\ 
and using the triangular inequality we obtain 



M 



M 



Sk 



MSk 



Since S'fc(r) = for \r\ < k, we can estimate this last term as follows: 



= / M 



Sk 



"Ipn'- 



where we recall that £k = {Q G ^ ; iV'nlQ)! ^ fc-M(Q)}. According to the Holder inequality, for 
all p > 1, denoting by q the conjugate of p (i.e. such that ^ + ^ = 1) and using the estimate (|4.7p . 
we obtain 



MS. 



"Km J 



M 



1/9 



< 



M 
M 



2p-^ 1/p 
2p\ 1/p 



We thus control the L^^-norm of MSk{^n/M) using his L^-norm . But this L^^-norms can itself 
be controlled, for an adapted value of p by the Hl^-noim. In fact, using the weighted Sobolev 
embedding (see Lemma I2.5P there exists p > 1 for which we have the inequality 



(I-IMI) 



2p\ 1/p 



< 



< 



Hi 



M*(t) 

MSk 





2 

+ 




1 


\m) 


A/,0 




J 2 



We deduce a control on the L|j-norm of MSk{ipn/M) using his HIj Q-norm: 



(ln(l + 



(ln(l + fc))2/'? 



< 



MSk 



\M 



HI 



Ll, ^ (ln(l + fc))2/9 

that is a control of the form WMSkitpn/M) ILa ^ -A{k)\iMSk{'4'n/M^ 11 „i where A{k) tends to 

when k tends to +oo. Hence, we obtain the following estimate for the term B of the left hand side 
of equation (|4.ip : 



H\ 



(4.9) \B\^\\^^U^^,,)(k + A{k)\MSk{^) 

• The last term of the equation ()4.ip is controlled as follow 

|C| = l(/,'^)l<llv'lk;, =^ 



MSk 



\M 



HI 



(4.10) 



M 



\m) 



< 



M 



H\ 



The previous estimates (j4.8p . (I4.9P and (|4.10p enable the deduction, from equation (|4.ip . for all 
e N, of the following inequality 



Since ^(fc) tends to when k tends to +oo, it possible to obtain for a sufficiently large fc, the 
inequality 

(4.11) ^^'^(^) 
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Estimate of MTk{rpn/M) in Jf^^Q-norm - Choose now (p = MTk{i)n/M) - M MTfc(V'„/M) 
as a test function in equation (|4.2p (according to Lemma 12.71 we have cp G H^jq). As for the 
estimate of MSk{ipn/M), we study each of three terms, named A, B and C as previously, present 
in equation (|4.2[) . 
• The first is written 



A 



M 



AM 



For the second term, we proceed as follow: 

\B\ ^ 

< Ik 



M 



Tn 



\M 



But for \ipn/M\ ^ k we have V(Tfe(V'„/M)) = whereas for \ipn/M\ < k we clearly have \ipn\ < kM 
and consequently, according to the Cauchy-Schwarz inequality we obtain 



\B\ < Ikl 



kM 



k^M, 



^k\\K\\L^^n)\\Mn(^^) 



M 



V Tfe 



M 



The last term is treated like those of the previous estimates: 



\c\ = \{f,^)\<\\Mn[^) 



Hi 



These three estimates give (note that this estimate depends on k, but that k has been fixed) 



(4.12) 



< 1. 



Estimate of ipn in ^ - Since for all fc G N we have St+Tk = id we obtain 



MSk 



m) 



Hi 



+ 













Using the estimates (|4.1ip and (|4.12p we deduce that for all n e N we have 

(4.13) < 1- 

Convergence of the sequence {i/'nlneN - According to the estimate (I4.13P the sequence 
{ipn}n<£N is bounded in ff q. According to the Lemma [2^ a subsequence of the sequence {'0n}neN 
(always denoted by {"i/'nlnGw) admits a limit ip, weak in q and strong in L\^. In order to perform 
the limit in equation (|4.2p . it is sufficient to prove that the sequence {MT„('0„/M)}„gN tends to V' 
in L\.j. We obtain 



MT, 



Li, 
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However the application T : K ^ M is 1-lipschitz and we have 

ii"-"(i?)-«-»(i)ii;^/.«i-»(i?)--"(i)r 

«/„"i^-sr=«*'"-«v 

which proves that ||MT„(-0„/Af) — AfT„('(/;/M)||^2^ tends to when n tends to +00. As regards the 
other term, the Lebesgue convergence dominated theorem directly affirms that \\MTn{ip/M) — ip\\i^2 
also tends to when n tends to +00. Finally, it was shown that the sequence {MTn{ilJn/M)}neN 
converges to ip in L^j and consequently that -0 is a solution of 



4.2. Uniqueness proof in Theorem 

Main steps for the uniqueness proof - To prove uniqueness, we proceed as follows: We start 
by introducing the dual problem. It is shown that this dual problem admits a solution by using the 
Schauder topological degree method. Then, by using the existence both problem and its dual, we 
deduce uniqueness from these two problems. 

Introduction of the dual problem - For g € ff^/ consider the elliptic partial differential 

equation 

(4.14) - div(A/v(^)) - Mk • v(^) = g on n 

and we look for a solution (f) G q to this equation. 

A compact application for the dual problem - For cj) e q we have Mk ■ V{(f)/M) E L\j C 
since \\Mk ■ V(0/Af)||^2^ ^ HKH^oofQ-) ||(/)|| jifi^ ^. Since the operator (p — div(Af V((/j/A/)) is 
coerciv in H\j q, there exists thus a unique solution </> = G(0) G H\.[ q such that for all f G H\j q 

This defines an application G : -ff If ~^ ^li o- quite easy to see that G is continuous; indeed, 
if 4>n tends to </> in i?|f q then MK-V{(l}n/M) tends to Mk-\/{(I)/M) in H~j^^ (more precisely in L\j). 
Thus div(MV((/)„/M)) tends to A\Y{MV{(j)/M)) which imphes that = G{^n) tends to = G{4)) 
in i?M,o- 

We will now prove that G is a compact operator. Suppose that the sequence {4>n\nen is bounded 
in Hj^ then {Mk ■ V((/)„/M)}„gN is bounded in H^^ so that, using ip = G((/)„) = as a test 
function in the equation satisfied by we get using the Lemma 12.101 



Mk-V{ 



(t?)lL-)"^""^M,0' 



which implies that the sequence {(pn}n£N is bounded in H^j q. Using the Lemma 12. 4[ up to a 
subsequence, we can thus suppose that {0„}„gN converges a.e. on fl and is bounded in Lj^. Let 
(n, m) G N^; subtract the equation satisfied by (f>m to the equation satisfied by (/>„ and use (p = 
i'n — 4'm as a test function, this gives using the Lemma l2. 101 again 



2 < 



M 



< ||0n - <P'm\\Ll 
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From the strong convergence of {(f)n}nefi to in L\j we deduce that the sequence {4>n}neti is a 
Cauchy sequence in H^j q and converges in this space. We deduce that the apphcation G is compact. 

Existence result for the dual problem using the Leray-Schauder topological degree - 

According to the Leray-Schauder topological theory (see the founder article of J. Leray and J. 
Schauder |16j ) since the operator G introduced with equation (j4.15p is a compact operator, to prove 
that it has a fixed point, we just have to find R > such that for all s S [0, 1] there exists no 
solution of — sG{(j)) — satisfying o ^ ^' 

Let s e [0, 1] and suppose that (f> e HIj g satisfies 4> = sG{(t>). We have for all Lp E H\,[ q 



(4.16) ( MV(^) . V(^) - . 1^ . V(|) = (sg, ^) 



Using the "non-dual" problem (see the existence proof of Theorcm l3. 21 where we obtain an existence 
solution of equation (|4.2p ). we know that for all / e H^^ there exists at least one solution ip £ H\,^ q 
such that for all if S H\j g 



Moreover, according to estimate (|4.13p there exists Ci € R+ such that for all / G ^a"/ with 
ll/ll^-i ^ 1 and for all s G [0, 1] we have HV'llffjj =^ Ci- We can verify that this constant C\ depends 
only on and can be selected independently on the function / when ||/||jLf-i ^ 1. In addition 

according to the estimates obtained in the existence proof of Theorem 13. 21 this constant Ci depends 
on IIskII/^oo^q) ^ ||'«||l°°(o) a-nd consequently the constant Ci can also be selected independently 
of s. 

By taking = in the equation (|4.17p satisfied by ijj and ip = ip in the equation (I4.16P satisfied 
by (f), we obtain 

(/,0) = (sg.ij) < s\\g\\^-iGi < Mh-^C'i C2. 

M M 

Since this inequality is satisfied for all / G H^l such that ||/||^-i ^ 1, we deduce that H?!'!!//!^ ^ ^ C2. 
Now take R — G2 + ^- We have just proven that, for any s E [0, 1], any solution of — sG{(t>) — 
satisfies ^' ^^^^ ^-^ Leray-Schauder topological degree theory, the application G has 

a fixed point, that is to say a solution of ()4.15p . 



Uniqueness - Since the equation (13. ip is linear, it is sufficient to prove that the only solution 
of ()3.ip without source term, i.e. taking / = 0, is the null function. Let ip he a, solution of ()3.ip 
with / = and let ^ be a solution of ()4.14p with g ~ sign(^/;) G . By putting = as a test 
function in the equation (|3.ip satisfied by ip and f = ip ds a, test function in the weak formulation 
of the equation ()4.14p satisfied by p, we respectively get 

We deduce that {sign{ip), ijj) — 0, that is to say /^^ \ip\ ~ and then tp = 0. 



Remark 4.1. — A similar reasoning gives the uniqueness of the solution of the dual prob- 
lem (1133 ■ 
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5. Application to fluid mechanics 

5.1. The FENE model for dilute polymers 

A natural framework where vectors fields strongly explode at the boundaries of a domain is the 
framework of the modeling of the spring whose extension is finite (that is physically realist). In fluid 
mechanics, such an approach is used to develop polymer models in solution. It is this point of view 
which we have chooses to present in order to illustrate the preceding theoretical study. 
The simplest micro-mechanical approach to model the polymer molecules in a dilute solution is the 
dumbbell model in which the polymers are represented by two beads connected by a spring. The 
configuration vector Q describes the orientation and the elongation of such a dumbbell lUj . The 
force of the spring is governed by some law that should be derived from physical arguments. We 
choose here the popular FENE model, in which the maximum extensibility of the dumbbell is fixed 
at some value determined by the dimensionless parameter £ and the spring force takes the simple 
form ^ 

^^^^ " T~|Q|W 

The configuration vector Q depends on time t and macroscopic position of the dumbbell x in the 
flow. Moreover, it satisfies the following stochastic differential equation (see [H] for details): 



(5.1) dQ=((Vu)^.Q-^E(Q))dt + ^^dW, 

where the 2-tensor (Vu)-^ is the transposed velocity gradient, I?e is a dimensionless number called 
the Deborah number (linked to the relaxation time of the fluid) and W is the Wiener random process 
that accounts for the Brownian forces acting on each bead. Equation (|5.ip should be understood as 
the Ito ordinary stochastic differential equations along the particle paths since the dumbbells'centers 
of mass are supposed on average to follow the particules of the solvent fluid. 

As is well known (see Section 3.3 of [21]), every Ito ordinary stochastic differential equation can be 
associated with a partial differential equation for the probability density function ip(t, x, Q) of the 
random process Q(t, x). In particular, equation (|5.ip implies the following, also called Fokker-Planck, 
equation for Lp{t, x, Q): 

(5.2) t + " • ^'^^ - ^^Q^ - dive (K(Vu)" ■ Q - ) . 

In certain modes the dominating terms correspond to the terms of the right-hand side member of 
equation (|5.2p . It is the case, for instance, when the flow is supposed to be thin, see [5]. In these 
configurations, the distribution ip can be seen like depending only on Q (to be rigorous, tp also 
depends on time t and on the macroscopic position x, via the presence of the gradient Vu(t, x) but 
these dependences can be seen as parameters) and the equation (j5.2p is approached by the following 
Fokker-Planck equation on f{Q): 

(5.3) - A</3+ div(<^F) = inB{0,£), 

with F = 22?e (Vu)-^ • Q — E. It exactly corresponds to those studied in the first part of this paper 
(see equation (jl.ip ) in the case fl = B{Q,£) and without source term: / = 0. 

Although it is wished that the solution (p{Q) cancels for values Q such that |Q| — £ (i.e. we 
wishe that the maximum length of the springs is £ and that there is no spring of this length), the 
classical framework of the Theorem 11.21 does not correspond to this equation provided with the 
homogeneous Dirichlct boundary conditions. In fact, the force F is not sufficiently regular: we have 
F ^ L'^{B{0,£)). Roughly speaking, the FENE model takes into account the finite extensibility of 
the polymer chain, through an important explosive force when |Q| tends to £. 
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_ - ■ ^ n-; 




Figure 5.1. On the left, we have drawn the physical domain of a real flow for a diluted polymers 
solution. From the microscopic point of view, the polymer chains are identified to independent 
mass-springs systems (called dumbbells). The orientation and the length of each dumbbell 
is governed by a quantity (denoted by ip in this paper) distributed in a ball whose radius 
corresponds to the maximum extension of the spring. On the right, the colors correspond to the 
various probabilities that dumbbell be in the given position. For instance, the drawn dumbbell 
is the dumbell which has the most chance to be present (with its "symmetrical" compared to 
the center of the ball) . 



On the other hand, this force F perfectly corresponds to the principal result shown in this article 
(see Theorem l3.3p . More precisely, the vector field Q G B{0, £) ^ 2'De (Vu)-'" • Q is clearly bounded 
and we can write the "explosive" term E as follow: 



E = -VV with 



nQ) = yin(i- 



IQP 



To make appear the maxwellian function M as it is used in this paper, we write 



(5.4) 



M(Q) 



1 



IQI 



5(0/) 



(l - \Rf/t 



dR 



From Theorem [331 we deduce that if £ > \/2 then for any p e M there exists a unique (weak) solution 
of the Fokker-Planck equation (|5.3p such that /g^Q ^-j f = P- According to H.C. Ottinger [3T], the 
number £ roughly measures the number of monomer units represented by a bead and it is generally 
larger than 10. The assumption £ > is not constraining from the physical point of view. In 
fact, according to H.C. Ottinger [21| . the number £ roughly measures the number of monomer units 
represented by a bead and it is generally larger than 10. 

Moreover, impose the quantity J^^q f physically corresponds to given the density of the polymer 
chains. Hence this condition is relevant for the studied problem. 

Remark 5.1. — If the tensor (Vxu)-'" is replaced by its anti-symmetric part ^(VxU — (Vxu)"^) 
in the force term F then we get the so-called co-rotational FENE model. This case corresponds to 
a particular cases presented page [51 f — M is a trivial solution of equation (|5.3p (see [51 [T5]V 



5.2. Numerical results 

In this subsection, we present numerical result for the Fokker-Planck equation (jl.ip for a confinement 
vector field F coupled with the normalization condition J^^ip = p, p € M, and then we apply the 
algorithm in the framework of fluid mechanics. The main difficulty to obtain a numerical scheme for 
the Fokker-Planck equation within the normalized condition is to treat this normalized condition 
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since the equation is not numerically difficult itself. Preci sly, this condition is implemented by 
penalization. For simulation, we use the FreeFemH — h progran ^ 
of the problem and finite elements method. 



(7) 



which is based on weak formulation 



In the fluid mechanics context, we want to observe the distribution of the orientation dumbells in 
a dilute polymer under shear (for instance with a given stationary velocity flow given of the form 
u(a;i,a;2) = (72:2,0), 7 G M, in the 2-dimensional case). For simplicity, we make the presentation with 
the 2-dimensional model. According to the previous subsection, the searched distribution satisfies 
the following Fokker-Planck equation 



(5.5) 



A(^+ div((^F) = onS(0,£), 



where the vector field F is given by 



e B{0,£) I — > 2Vej 



Q2 




IQI 



Qi 
Q2 



Moreover, the solution must be satisfy the relation J^^^ ^^tp ~ p. Notice that if we have a solution 



such that /g(o 'Z' = 1 then, by linearity, the function Lp — pip is a. solution such that 



B(0,- 



In the numerical test, we always take p — 1. The only two parameters which are interest are the 
product Ve 7 and the coefficient £ which corresponds to the maximal elongation of the dumbells. 
Without shear (that is for 7 = 0), a trivial solution of the Fokker-Planck equation (|5.5p exists: 
it is the maxwcllian M (see its expression (|5.4p ). For three characteristic maximal lenghts of the 
dumbells {£ — 2, £ — 5 and £ — 10), we have been represented this maxwellian on the figure [521 






Figure 5.2. Solution without shear for different maximal lengths of dumbells: £ = 2, £ = 5 and 
£^ 10. 

To observe the influence of the shear on the distribution, taking De = 10, ^ = 5 and different values 
of the shear coefficient 7 £ {0.1; 0.2; 0.5; 1}. The for results are descibe on figure [5?3l 



Figure 5.3. Shear infiuence on the distribution of the length of dumbells. The more the shear 
is raised, the more the dumbbells tend to lengthen in the direction of the fiow. The four figures 
above correspond (from left to right) to the values 7 = 0.1, 7 = 0.2, 7 = 0.5 and 7=1. 



see |http://www.freefem.org/ff[ f-|- 
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